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Abstract 

We searched for a sidereal modulation in the rate of neutrinos produced by 
the NuMI beam and observed by the MINOS far detector. The detection of 
such harmonic signals could be a signature of neutrino-antineutrino mixing 
due to Lorentz and CPT violation as described by the Standard-Model Ex- 
tension framework. We found no evidence for these sidereal signals and we 
placed limits on the coefficients in this theory describing the effect. This is 
the first report of limits on these neutrino-antineutrino mixing coefficients. 
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1. Introduction 

The sensitive interferometric nature of neutrino oscillations offer interest- 
ing and novel ways to search for Lorentz invariance violation (LV) and CPT 
violation (CPTV) that may arise from physics at the Planck scale, 10 19 GeV. 
The SME is effective field theory that contains the Standard Model, General 
Relativity, and all possible effective operators for Lorentz violation |lj. In 
the SME, LV and CPTV could manifest themselves at observable energies 
through a dependence of the neutrino oscillation probability on the direction 
of neutrino propagation with respect to the Sun-centered inertial frame. An 
experiment that has both its neutrino beam and detector fixed on the Earth's 
surface could then observe a sidereal modulation in the number of neutrinos 
detected from the beam. 

MINOS is such an experiment [2]. It uses Fermilab's NuMI neutrino 
beam [3] and two neutrino detectors - a near detector (ND) that is 1.04 km 
from the beam target and a far detector (FD) that is 735 km from the beam 
target. Both detectors are magnetized to approximately 1.4 T, allowing 
for the identification of fi~(fi + ) produced in charged-current (CC) neutrino 
(antineutrino) interactions. In a series of recent papers jH [5], |6], MINOS 
searched for LV and CPT violation with neutrino and antineutrino data 
recorded in its detectors 0[B]. These analyses focused on detecting oscillation 
signals in neutrino and antineutrino flavor-changing reactions that vary with 
a sidereal period, as predicted by the Standard Model Extension (SME). 
MINOS found no evidence for LV and CPTV in these analyses. In the present 
paper, we extend the search for LV and CPTV using the timestamps of the 
neutrino data made available by the MINOS collaboration to look for sidereal 
signals that would indicate neutrinos changing to antineutrinos. Although 
there have been many searches for LV and CPTV with the SME [9], this is 
the first to probe the possibility of neutrino-antineutrino mixing. 

Denote the usual neutrino survival probability in the two-flavor approx- 
imation as Pvj}-*Vn ~ 1 — sin 2 (26> 23 )sin 2 (1.27Am 2 2 L/£), where 8 23 is the 
angle describing mixing between the second and third mass states and Am^ 
is the difference in the squares of those mass states. The energy of the neu- 
trino is E and the distance it travels is L. Then LV and CPTV that would 
cause mixing between neutrinos and antineutrinos introduces an additional 
perturbation term to the survival probability [10J, 

p = p(0) , p(2) (1 \ 
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where the perturbation term PvJ^v^ can be written 

P %1»» = L 2 ^P c + P Bs sm2u (B T (B + P Bc cos2u (B T (B (2) 

+ Pjr s sin 4o;©Te + P Tc cos 4w e T ffi j . 

Here L = 735 km is the distance from neutrino production in the NuMI beam 
to the MINOS FD [llj, tu e = 27r/(23 h 56 m 04.0982 s ) is the Earth's sidereal 
frequency, and T e is the local sidereal arrival time of the neutrino event. The 
parameters Pc, Ps a , Pb c , Pj^ s , and Pjr c contain the LV and CPTV information 
on neutrino-antineutrino mixing. They depend on the SME coefficients H°t 
and g°f, the neutrino energy, and the direction of the neutrino propagation 
in a coordinate system fixed on the rotating Earth [10J. As is clear from 
Eq.([2]), we need only the even harmonics 2o; ffi T ffi and 4o; e T e for this analysis. 
Further, since we are testing the theory with a harmonic analysis, there is 
no sensitivity to Pc- 

In the minimal SME, neutrino-antineutrino mixing v a -H- Uj, is controlled 
by the complex coefficients and g ^ (a = e, /i, r and b — e, ft, f). Sidereal 
modulations detectable by MINOS are produced by 6 of the H°r coefficients 
and 60 of the g"? coefficients [10]. For the coefficients, only three in- 
dependent pairs of flavor subscripts (ab = e/2, ef , /if) and two spatial com- 
ponents (a = X, Y) produce sidereal variations. For the g ^ coefficients, 
MINOS can access six pairs of flavor indices (ab = ee, eft, ef, \ift, \if, tT ) and 
10 different spacetime components (a(3 = XT, YT, XZ, YZ, ZX, ZY, XX, 
YY, XY, YX). 

For the neutrino propagation direction, we specify the direction vectors in 
a right-handed coordinate system defined by the colatitude of the NuMI beam 
line x = (90°- latitude) = 42.17973347°, the beam zenith angle 9 = 86.7255°, 
measured from the z-axis and which points toward the local zenith, and the 
beam azimuthal angle, <p = 203.909°, measured counterclockwise from the 
x-axis, which points south. 

2. Data Analysis 

This analysis uses a data set of neutrino interactions acquired from May, 
2005 through April, 2012 [7J [S] when the beam was running in a mode that 
produces 91.7% v^, 7% V^, and 1.3% v e + V e [12]. The interactions were 
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selected using standard MINOS criteria for beam and data quality. In addi- 
tion, the events were required to interact within the 4.2 kiloton FD fiducial 
volume. This selection enables MINOS to establish each event as a CC 
interaction by identifying the outgoing yT . The method is described more 
fully in [TJ E]. As in [U El E] , we focused on CC events to maximize the 
disappearance signal. There are a total of 2,463 CC events in this analysis. 

The analysis proceeded by first tagging each neutrino interaction with the 
local sidereal time (LST) of its spill - the GPS time of the spill, accurate to 
200 ns [13], converted to sidereal time. The event's LST was then converted 
to local sidereal phase, LSP = LSTx (w ffi /27r), a parameter in the range 
0-1. Event times were not corrected for their time within a 10 /is spill, 
an approximation that introduces no significant systematic error into the 
analysis. The data for each spill were binned into two histograms. In one 
histogram, we added any interaction from each spill; in the second histrogram 
we added the number of protons impinging on the target (POT) for each 
spill. These two histograms have 32 bins in LSP because the Fast Fourier 
Transform (FFT) algorithm we use to search for sidereal variations [H] works 
most efficiently for 2^ bins, where J\f — 1 is the number of harmonics to be 
computed in addition to a constant term. Since the highest harmonic in 
Eq. ^ is 4a; e T e , we adopted J\f = 5 as the binning. Each phase bin spans 
0.031 in LSP or 45 minutes in sidereal time. 

After adding all the events and POT in the data set into the two his- 
tograms, we divide them to get the rate of interactions per POT, [y^ CC- 
events) / (POT) as a function of LSP. This histogram is the one used to search 
for sidereal modulations in the neutrino rate. Figure [T] shows this histogram 
for our analysis. The statistically significant fit to a constant rate implies 
there are no sidereal modulations in the data sample. 

We next performed an FFT on the rate histogram in Fig. [T] and com- 
puted the power in the harmonic terms 2o; e T e and 4w ffi T ffi appearing in the 
oscillation probability, Eq.Q. Let S2 be the power returned by the FFT for 
the second harmonic term sin (2o; e T e ) and C2 be the power returned for the 
second harmonic term cos(2o; e T e ); similarly define S4 and C 4 . Then the 
statistics we used in our search are 



We added the powers in quadrature to eliminate the effect of the arbitrary 
choice of a zero point in phase at h LST. The results of the FFT analysis 
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Figure 1: The phase diagram of the CC neutrino event rate for the FD data. The mean 
rate of 2.31 ± 0.05 events per 10 15 POT is superposed. For the fit, x 2 /ndf = 23.7/31. 

are given in Table [TJ Although only the even harmonics p-i and p^ are needed 
for this analysis, we include the power in the odd harmonics p± and p% as a 
check on whether alternative theories on LV that put power into these modes 
are realizable. 

Table 1: Results for the Pi, - ■ ■ ,Pi statistics for the histogram shown in Fig. [T] The third 
column gives the probability, Vp, that the measured power is due to a noise fluctuation. 



Statistic 


p(FFT) 


v T 


Pi 


0.928 


0.65 


P2 


0.574 


0.89 


P3 


1.388 


0.48 


P4 


1.223 


0.53 



As described in [5], we determined the statistical significance of the har- 
monic powers p(FFT) in Table [I] by using the data themselves to construct 
10 4 simulated experiments without a sidereal signal. To construct a simu- 
lated experiment we used the sidereal time distribution of the beam spills to 
determine a new, randomly selected LSP for each interaction in the data. Us- 
ing this randomly chosen phase, we added each interaction to one histogram. 
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Similarly, for each spill we selected a random LSP and added the number 
of POT in the spill to a second histogram. The FD records fewer than one 
interaction per spill, so the spill and interaction LSP for these simulated 
experiments can be selected independently without introducing any biases 
in the procedure. As for the real neutrino data set, the histograms had 32 
bins in LSP. We then divided the events histogram by the POT histogram 
to obtain the rate histogram for the simulated experiment. We repeated 
this process 10 4 times. This procedure ensures that the spill times are dis- 
tributed properly in LSP. Most importantly, it scrambles any coherent signal, 
if present, in the data and eliminates it in the simulated experiments. 

We next performed an FFT on each simulated histogram of (events/POT) 
and computed the power in the four harmonic terms (a>®T®, . . . ,4c%T©). 
Fig. [2] shows the distribution of p\, . . . ,p± for the 10 4 simulated experiments. 
The distributions for pi, . . . ,p^ are quite similar over the full range of ob- 
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Figure 2: The distributions for the quadratic sum of powers Pi, - ■ ■ ,Pa from the FFT 
analysis of 10 4 simulated experiments without a sidereal signal. The inset shows the 
distribution for P2 with a fit to a Rayleigh distribution with a = 1.49 superposed. 

served powers. As seen in the inset to Fig. |2j the P2 distribution is well 
described by a Rayleigh distribution, as expected in these simulated experi- 
ments since they were constructed so that the power in the sine and cosine 
terms are uncorrelated and normally distributed. 

We assume that the threshold for signal detection in any harmonic is the 
quadratic power p(FFT) that is greater than 99.7% of the entries in its respec- 
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tive Pi, ■ ■ ■ ,Pi histogram. The values of p 9 g.7%(FFT) for all four harmonics 
lie between 3.40 - 3.50, with p 2 = p 4 = 3.45. We adopt p 99 . 7% (FFT) = 3.45 
as the 99.7% confidence level (C.L.) for the probability that a measured 
quadratic sum of powers for any harmonic was not drawn from a distribution 
having a sidereal signal. As neither p 2 nor p 4 exceed our detection threshold, 
we conclude that there is no evidence for a sidereal modulation resulting from 
neutrino- antineutrino mixing as predicted by Eq.([2]) in the neutrino data set. 
This result is consistent with the conclusion that the constant rate fit to the 
data in Fig. [T] implies there is no sidereal modulation in the data. 

We determined the minimum detectable sidereal modulation for this anal- 
ysis by injecting a sidereal signal of the form A sin(u; e T), where A is a fraction 
of the mean event rate, into a new set of 10 4 simulated experiments and re- 
peating the FFT analysis. We found that every experiment gave p > 3.45, 
or exceeded the detection threshold, when A = 13% of the mean rate. 

In Table [T] we also give TV, the probability that the power measured in a 
particular harmonic is due to a noise fluctuation. In other words, Vjr is the 
probability of drawing a value at least as large as found for pi, . . . ,p 4 from 
the parent distribution in Fig. [2] of the powers for numerical experiments 
without a sidereal signal. The large values for Vj also make the case that 
the powers we measure are the result of statistical fluctuations. 

We investigated the sensitivity of our results to several sources of system- 
atic uncertainties. In the previous MINOS analyses [U[5], the degradation of 
the NuMI target was observed to cause a drop in the event rate. Further, the 
NuMI target was replaced multiple times during the data taking period of 
this analysis which could also have affected the event rate. We investigated 
how these changes in the event rate would affect the detection thresholds and 
found this source of systematic uncertainty to be negligible. 

Day-night effects are another potential source of systematic uncertainty 
for a data set like ours that does not have uniform coverage throughout a solar 
year. We find the mean rate during the day to be 2.32±0.06 events/10 18 POT 
and the mean rate at night to be 2.21 ±0.07 events/10 18 POT. Since there is 
no statistically significant difference in these rates, we conclude that diurnal 
effects cannot be masking a true sidereal signal in the data. 

There is a known uncertainty of ±1% in the recorded number of POT per 
spill [15J . We verified this uncertainty could not introduce a modulation that 
would mask a sidereal signal as done in [6] . We further determined that long 
term drifts in the calibration of the POT recording toroids of the size ±5% 
over six months would not alter the detection threshold. We thus conclude 
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that POT counting uncertainties are not masking a sidereal signal. 

Compton and Getting [16j first pointed out that atmospheric effects can 
mimic a sidereal modulation in the event rate if there were a solar diurnal 
modulation that beats with a yearly modulation. As in [6], we found the 
amplitude of the potential faux sidereal modulation would be only 1.6% of 
our minimum detectable modulation and therefore would not mask a sidereal 
signal in these data. 

3. Limits 

We determined the confidence limit for a particular SME coefficient by 
simulating a set of experiments in which we inject an LV signal into a Monte 
Carlo simulation that does not include a sidereal modulation. In this simu- 
lation, neutrinos are generated by modeling the NuMI beam line, including 
hadron production by the 120 GeV/c protons striking the target and the 
propagation of the hadrons through the focusing elements and decay pipe 
to the beam absorber. It then calculates the probability that any neutrinos 
generated traverse the FD. The FD neutrino event simulation takes the neu- 
trinos from the NuMI simulation, along with weights determined by decay 
kinematics, and uses this information as input into the simulation of the 
interactions in the FD. We inject a sidereal signal in the simulation by calcu- 
lating the survival probability for each simulated neutrino based on Eq. ([TJ 
using a chosen value for the magnitude of the SME coefficient, the energy of 
the simulated neutrino, and the distance the neutrino travels to the FD. We 
start by setting the SME coefficient to zero and then continually increasing 
its magnitude until either p2 or crosses the detection threshold of 3.45. 
We took this value of the SME coefficient to be its 99.7% CL. upper limit. 
The SME coefficients can be either positive or negative, so this limit is the 
lower limit on the SME coefficient in the case it is negative. We repeated 
the simulation for this coefficient 250 times and averaged the resulting value 
of the coefficient. Each of the simulated experiments contained the same 
number of interactions as the data set. 

The 99.7% CL. limits on the SME coefficients are given in Tableland 
Table |3j By setting all but one SME coefficient to zero to determine its 
confidence limit, our method is based on the premise that our null detection 
does not result from fortuitous cancellations of SME coefficients that hide a 
signal of oscillation terms in Eq. ([2]). Since the number of SME coefficients 
is large, this could be an issue. This issue was explicitly addressed in [B] and 
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Table 2: Upper limits to the modulus of the IZe parts of the coefficients in [GeV] 
describing ncutrino-antincutrino mixing in the SME theory. The limits were determined 
from the 99.7% threshold for a positive sidereal signal in the neutrino rates in the MINOS 
FD. The IZe and 1m components of these SME coefficients are computed identically and 
are equal. 



a 




\Ke{H«,)\ 


\Ke{H^)\ 


X 
Y 


3.3 x 10" 21 
3.3 x icr 21 


3.4 x 10" 21 
3.6 x icr 21 


i.i x icr 22 
l.o x icr 22 



such a scenario was shown to be quite improbable for the flavor changing co- 
efficients investigated in those previous analyses. We assume that conclusion 
holds for this analysis also. 

4. Summary 

We have presented a search for the Lorentz and CPT violating sidereal 
modulations in the observed neutrino rate predicted by the SME theory 
for neutrino-antineutrino mixing in the MINOS far detector. We found no 
significant evidence for this predicted signal. When framed in the SME the- 
ory [T71 [IS] , this result leads to the conclusion that we have detected no 
evidence for Lorentz invariance violation, a result consistent with the anal- 
yses in [U El E]- We computed upper limits for the 66 SME coefficients 
appropriate to this analysis. These limits are the first to be calculated for 
the SME coefficients governing neutrino-antineutrino mixing. 
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Table 3: Upper limits to the modulus of the TZe parts of the coefficients describing 
neutrino-antineutrino mixing in the SME theory. The limits were determined from the 
99.7% threshold for a positive sidereal signal in the neutrino rates in the MINOS FD. The 
TZe and Tm components of these SME coefficients are computed identically and are equal. 



a(3 


\ne(~g e f)\ 


\Ke(g e f)\ 


l^eaf)| 


\i**<3$)\ 


\Ke{g$)\ 




XT 


7.6 x 10- 


22 


7.6 x 10- 


22 


8.2 x 10- 


22 


8.9 x 10" 


-24 


8.6 x 10- 


24 


1.8 x 10" 


22 


YT 


7.6 x 10- 


22 


7.6 x 10- 


22 


7.6 x 10- 


22 


8.6 x 10" 


-24 


8.4 x 10- 


24 


1.8 x 10- 


22 


XZ 


1.2 x 10- 


21 


1.2 x 10- 


21 


1.2 x 10- 


21 


1.3 x 10- 


-23 


1.3 x 10- 


23 


2.9 x 10- 


22 


YZ 


1.2 x 10- 


21 


1.2 x 10- 


21 


1.2 x 10- 


21 


1.3 x 10- 


-23 


1.4 x 10- 


-23 


2.9 x 10- 


-22 


ZX 


1.0 x 10- 


21 


1.0 x 10- 


21 


1.0 x 10- 


21 


1.2 x 10- 


-23 


1.1 x 10- 


-23 


2.4 x 10- 


-22 


ZY 


1.0 x 10- 


21 


1.0 x 10- 


21 


1.0 x 10- 


21 


1.1 x 10- 


-23 


1.2 x 10- 


-23 


2.4 x 10- 


-22 


XX 


2.0 x 10- 


21 


2.0 x 10- 


21 


2.0 x 10- 


21 


2.3 x 10- 


-23 


2.3 x 10- 


-23 


4.8 x 10- 


-22 


YY 


2.0 x 10- 


21 


2.1 x 10- 


-21 


2.1 x 10- 


21 


2.2 x 10- 


-23 


2.2 x 10- 


-23 


4.8 x 10- 


-22 


XY 


2.0 x 10- 


21 


2.0 x 10- 


21 


2.1 x 10- 


21 


2.3 x 10- 


-23 


2.3 x 10- 


23 


4.8 x 10" 


22 


YX 


2.0 x 10- 


21 


2.0 x 10- 


21 


2.0 x 10- 


21 


2.2 x 10- 


-23 


2.2 x 10- 


23 


4.8 x 10- 


22 
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